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Despite the intensive efforts for determining the mechanism that causes high-temperature super- 
conductivity in copper oxide materials, no consensus on the pairing mechanism has been reached. 
Recent advances in high resolution angle-resolved photoemission spectroscopies have suggested that 
a sizeable electron-phonon coupling exists as the principal cause for kinks in the dispersion rela- 
tions (energy versus wave vector) of the electronic states. Here, we report on a systematic study 
about the influence of the electron-phonon coupling parameter "A" in the electronic quasiparticle 
^^ ' dispersions along the nodal direction for La2-xSrxCu04, Bi2Sr2CaCu20s+x and Y Ba2Cu30fi+x- 

yJ , Information about the dressing of the charge carriers, i.e., on the enhancement of the effective 

r^^ ■ mass and the strength of the coupling mode, is obtained as a function of the doping concentra- 

r~^, ' tion, temperature, momentum and energy from the kink dispersion in the (0-0)-(7r,7r) direction 

of momentum-space avoiding the complications of the d-wave superconducting gap. Our analysis 
shows a remarkable agreement between theory and experiment for different samples and at different 
doping levels. This includes our recently introduced theoretical model to adjust the experimental 
H I data of the fermionic band dispersion, emphasising the necessary distinction between the general 

electron mass-enhancement parameter A* and the conventional electron-phonon coupling parameter 
A. In LSCO, the coupling constant A, calculated consistently with the nodal kink dispersions, repro- 
5-H ' duces the observed critical temperatures Tc, the gap ratio 2Ao/kBTc, and other parameters which 

r^, have been studied from several equations. It will be concluded that the strong renormalisation of 

^ . the band structure can be explained in terms of the phonon coupling mode, and must therefore 

be included in any microscopic theory of superconductivity, even for those materials in which the 
contribution to the pair formation can be less dominant. Nevertheless, it seems unavoidable to 
consider additional mechanisms that justify the higher critical temperatures observed in BSCCO 
and YBCO samples. 
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CN ■ I. INTRODUCTION 

> 

t — ' One of the most relevant milestones in solid state physics has been undoubtedly the discovery of the high super- 

r~^ . conducting transition temperature (Tc) in LaBaCuO ceramics, prompting an intense activity in the field of synthesis 
^^ ' and characterisation of copper oxide materials (cuprates) [l|. The chemical composition of cuprates is well charac- 
ly-T , terised by a layered crystal structure with one or more Cu02 planes per unit cell, which are responsible for the low 
f— ^ ' lying electronic structure, and also presumably, the key for understanding the quasiparticle properties involved in 
f^ ] the mechanism of superconductivity, i.e.: the origin of the coupling mode which binds two electrons (holes) in the 
T— I . formation of Cooper pairs. Nevertheless, despite over two decades of experimental and theoretical efforts to elucidate 
the pairing mechanism and to understand the many-body properties influencing the electron dynamics in cupratc su- 
perconductors, there is no consensus on the strength of the electron correlations and the nature of the boson coupling 
modes. Nowadays, with the advent of improvement in instrumental resolution, sample quality, as well as theoretical 
development, the relevance of electron-boson coupling mechanisms on the electronic dynamics of cuprates is more 
and more accepted, because the many body effects arc supposed to be particularly strong (2|. In this paper, wc will 
review recent progress in the understanding of electrons interacting with bosonic modes, such as phonons, and discuss 
the results of recent high-resolution angle-resolved photoemission spectra (ARPES) on cuprate superconductors with 
the aim of studying its possible connection with the occurrence of high-Tc. In order to understand the relevance of 
the ARPES technique, we recall that the photoemission process results in both an excited photoelectron and a pho- 
tohole in the final state. Related to this observation, it becomes a useful probe of the related scattering mechanisms 
contributing to the electrical transport in different materials. Unlike other probes of the transport properties, the 
ARPES technique has the advantage of momentum resolving. Along this line, we want to note that the single-particle 
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scattering rate measured in ARPES is not identical to the scattering rate measured in transport studies themselves. 
Nonetheless, direct proportionality between them has been established 0, |3|- 

In conventional metals, the phonon coupling mode has long been recognised as the main mechanism involved in the 
superconducting properties, and the strength of this interaction essentially determines the value of Tc. On the other 
hand, in the high-temperature superconductors (HTSC) the experimentally determined d— wave pairing introduces 
considerable complications in the theory even when other coupling modes are considered [5|. Fortunately, with the 
appearance of a new era of analysers with improved resolution both in energy and momentum, in angle-resolved 
photoemission spectroscopy (ARPES) as well as in inelastic neutron scattering (INS) and x-ray scattering (IXS) 
experiments, the controversy on the influence of the anisotropic character of the superconducting gap in the electron 
properties can be avoided, by analysing preferential directions within the Cu02 planes or nodal directions. This 
provides a smart solution if one is merely interested in identifying the energy modes [6|. In fact, one of the most telling 
manifestations of the electron-phonon coupling is a mass renormalisation of the electronic dispersion at the energy 
scale associated with the phonons. This renormalisation effect is directly observable in the ARPES measurements as 
a low-energy excitation band in the dispersion curves of photoemitted electrons, known as kink 0, |8|. This feature, 
so far universal in HTSC, has been regarded as a signature of the strength of the boson mechanism which causes the 
pair formation in the superconducting state. In fact, all the interactions of the electrons which are responsible for the 
unusual normal and superconducting properties of cuprates are believed to be represented in this anomaly |9| . This 
has prompted an intense debate about the nature of the coupling mode involved in the density of low-energy electronic 
excited states in the momentum-energy space of cuprates, and its influence on the emergence of the superconducting 
state [i-il. 

So far, the origin of high-temperature superconductivity and the nature of the bosons involved remains controversial 
mainly because experiments can only be used to determine an approximate energy of the mode and this energy is 
close to both the optical phonons [3, [3 [ill, Il9l - [2l| and magnetic excitations [l^, llJ, l22l. |25| - |27| . Related to this, the 
energy scales of the optical phonons in the CUO2 planes are similar for electron and hole doped cuprates, while the 
magnetic mode for electron doped HTSC is found to be much smaller [2^[23|. Along this line, recently S. R. Park et 
al. j23| have demonstrated that the magnetic resonance mode can not explain the ARPES spectra in several electron 
doped HTSC systems showing a clear support for the electron-phonon coupling. In addition, J. Graf et al. |20| have 
recently reported the first evidence of an anomalous dispersion of the Cu-0 bond stretching phonon mode in a Bi 
cuprate supporting the idea that strong electron energy dispersion measured by ARPES corresponds to the Cu-0 
bond stretching phonon mode. On the other hand, considering that the magnetic resonance has not been detected in 
the single layer Bi2201 [30], magnetic modes should be ruled out as a general mechanism within a controllable theory 
of strong correlations. 

In an effort to clarify the influence of the phonon coupling mode (either weak or strong), in this work we are going 
to analyse the influence of the so-called electron-phonon coupling parameter on the electronic dispersion relations in 
several cuprate compounds. 

The paper is organised as follows. In Sec. [Hi we put forward some details about the conventional electron-phonon 
coupling theory focusing on the determination of the complex electron self-energy, and evaluating the influence of the 
electron-phonon spectral density calculated by different methods available in the literature on the superconducting 
properties of some cuprates (e.g., LSCO, YBCO, BSCCO). We concentrate on the use of approximate formulas to 
determine the superconducting transition temperature Tc, the ratio gap 2Ao/fcBTc, and the zero temperature gap Aq. 
The physical interpretation of the underlying approximations is also focused on. In Sec. IIIII we present a systematic 
study on the influence of the electron-phonon coupling parameter in the electron quasiparticle dispersion relations 
unravelling possible new physics. Comparison with the electronic dispersion curves measured by ARPES at the nodal 
point will be emphasised for the three families of cuprates. An empirical equation is incorporated for determining A 
from the doping level in LSCO which may be of interest for other works. Finally, Sec lIVl is devoted to discussing our 
results. The relevance of the electron-phonon coupling mechanism for the interpretation of the electron dynamics in 
HTSC will be concluded. Comparison with recent published material will be emphasised. 

II. ELEMENTS OF ELECTRON-PHONON COUPLING THEORY 

The purpose of this section is a brief reexamination of the conventional theory for electron-phonon coupling in 
condensed matter, with the aim of obtaining an appropriate interpretation of the so-called spectral function. This 
will envisage to extract some general consequences and then, check them from the direct interpolation between theory 
and experiment. 

Recall that, nowadays the simplest interpretation of microscopic superconductivity remains to be the Bardeen, 
Cooper, and Schrieffer's [31[ theory (BCS) which is based on the condensation of pairs of electrons (holes) into a spin 
singlet induced by an isotropic electron-phonon coupling. The formation of pairs gives rise to an excitation energy 



gap A in the electron density of states across the Fermi level {Ep) and its magnitude reflects the strength of the 
pairing interaction. In BCS the pairing function A has an isotropic s-wave symmetry, and the actual reduced gap 
value 2Ao/kBTc is used to classify a system as a weak- or strong coupling superconductor by comparison to the 
mean-field value 3.51. A significant accumulation of microscopic information on the fine structure found in tunnelling 
experiments at voltages above the forbidden energy gap was the most direct proof of the electron-phonon interaction 
as the coupling mechanism for superconductivity in conventional metals. However, along with the confirmation of the 
BCS predictions, the experimental results also indicated deviations from the weak-coupling BCS theory through the 
detection of large gaps as compared to the mean-field value. For the new generation of HTSC the detailed structure 
in the measured quasiparticle density of states, also revealed the need of one controllable strong-coupling theory. 
Moreover, many of the new families of superconductors discovered after the appearance of the cuprate compounds 
also exhibit highly anisotropic gaps, and transport and magnetic properties that cannot be directly explained by 
the simple BCS theory [32|- At present, there is no consensus on the microscopical theory involved in the HTSC 
materials and the influence of the pairing mechanism on the superconducting gap and the higher Tc- Nevertheless, 
assuming that the quasiparticle dynamics is governed by some electron-boson coupling which predominates both in the 
normal state properties and in the superconducting state, a wide number of properties can be satisfactorily explained 
[5l.l33j. Thus, having made these general considerations on the many-body problem underlying in the electron-phonon 
coupling theory, we are going to expose some technical aspects to be considered in the analysis of ARPES data and 
their relation with the superconducting properties. 

A. The Migdal-Eliashberg approach: spectral functions 

In the diagrammatic language, the above mentioned physical properties are formulated within the framework of the 
Fermi-liquid theory, where electron-like quasiparticles populate bands in energy-momentum space up to the cut-off at 
the Fermi energy. Lattice vibrations couple to electrons because the displacements of atoms from their equilibrium 
positions alter the band dispersions, either lowering or raising the total electron quasiparticle energy. 

Electron correlations are responsible of the formation of quasiparticles which are well deflned near the Fermi level. 
There, the so called vertex corrections can be a priori neglected because these can be shown to be reduced by the ratio 
between the phonon frequency (0 — 100 meV) and Ep (~ 1 — 10 eV). In other words, the presence of strong electron 
correlations mediated by the electron-phonon interaction is avoided and the multi-phonon excitations are reduced to 
the single-loop approximation or Migdal-Eliashberg (ME) approach (Fig.[T]d). In ordinary metals, this simple picture 
for the phonon-mcdiated interaction between electrons, or electron-phonon coupling, has long been known to be the 
pairing mechanism responsible for superconductivity. On the other hand, in the HTSC it has been suggested that 
other collective excitation modes mediate the pairing since the superconducting transition temperatures T^ are much 
higher than those of conventional superconductors. Nonetheless we want to emphasise that the feasibility of the ME 
approach in the HTSC is an issue still open. In fact, the recent observation of similar renormalisation effects in the 
energy-momentum space of HTSC has raised the hope that the mechanism of hi gh-r e superconductivity may finally 
be resolved, thus reviving the interest on the diagrammatic approach 0, Ht^j [2Jj ^3, [SJ] ■ 

Let us go into some more detail. The thermodynamic properties of the solid can be evaluated from the thermal 
Green's-function formalism where it is assumed that the electron-phonon interaction can be introduced by the relation 

G^^{k,iujn) = GQ^{k,iujn) - S(A:,ia;„) (1) 

with Gq related to the bare electron energy, i.e., the non-interacting electron Green's function {Go{k, ia;„) = (ia;„ — 
£k)~^) and iun standing for the so-called "imaginary Matsubara frequencies" [35|. 

Technically, the lowest order Fcynman graph approximation for the electron phonon interaction can be defined 
on a set of orthonormal functions (harmonic Fermi surfaces) that allow to consider the complexities of the d— band 
electron structure and phonons from the "angular" and "energy" components of the phase space [3g|, i.e.: fc-space is 
represented in terms of its harmonic components (J, e). Thus, in terms of this set 

xGj,{e\iuJn~i^,^) (2) 

where iV(0) — J^k^i^k) is the electronic density states on the Fermi surface, and the electron-phonon spectral 
density a'^F{J, J' ,e,e',v) represents a measure of the effectiveness of the phonons of frequency v in the scattering 
electrons from k(J,e) to k'{J',e'). Up to this point, only the harmonic approximation to the phonon propagator 
D{k — k',v) = 2u/{uj1 + v'^) has been considered. Nevertheless, we recall that the actual form of Eq. ^ is cumbersome. 



and requires of an aecurate determination of the electron-phonon spectral density from theoretical ab-initio calculations 
of the band structure of the crystal. Then, the ME approach is invoked. In order, to simplify Eq. © one assumes 
that it is possible to neglect the dependence on the energy surfaces (e,e') of the N{e')a^F{J,J',e,e',i') function. 
This allows us to omit the processes, violating the Born-Oppenheimer adiabatic theorem contained within the high 
order graphs. This is the ME "approach". Motivations for the use of this approach have been broadly discussed in 
the past (see Refs. [33| & (S^l)- 

So far, the self-energy has been described in the orthonormal space of "Fermi surface harmonics" which is to be 
transformed in the phase-space under consideration for k-spacc rotations according to their irreducible representations 
of the point group of crystal. Formally, this means that the spectral function a^F{JJ', v) is to be diagonal. It must be 
emphasised that this diagonal representation for the normal state will continue to hold in the superconducting state 
for the isotropic Cooper pairing or s-wave gap. However, other pairing schemes which break rotation symmetry are, in 
principle, possible. Thus, taking advantage that the ARPES measurements at the nodal direction, are not influenced 
by the anisotropy of the superconducting gap, we will refer to a (nondirectional) isotropic quasiparticle spectral 
density, defined as the double average over the Fermi surface of the electron-phonon spectral density a^i^(k,k', z^); 
i.e., 

«'^M = ]^ E I aL' P '5(^ - i^U,)Sie,)Sisk') (3) 

where, g^^., = [h/2Mi'^,-^^]^/^{'k\e-^,-^ ■ Vy|k') defines the matrix elements of electron-phonon interaction for electron 
scattering from k to k' with a phonon of frequency i^k„[(.' i^ ^ branch index). M stands for the ion mass, V is the 
crystal potential, e^,^ is the polarisation vector, and iV(0) — '^f^S{ek) represents the single-spin electronic density 
of states at the Fermi surface. As usual d{x) denotes the Dirac's delta function evaluated at x' = 0. In addition, 
note that ISk^'l^ ^^ inversely proportional to the number of charge carriers contributed by each atom of the crystal 
to the bosonic coupling mode. Therefore, an increase in the doping level, which causes an increment in the hole 
concentration of the Cu02 plane must be reflected in the coupling parameters as we will see in section Hill Moreover, 
recalling the outstanding feature of the theory of metals, that |3k]5./p vanishes linearly with |k — k'| when |k — k'| <C kp 
p38J , one would expect a linear disappearance of the coupling effect that gives rise to the nodal kink in the vicinity 
of the Fermi surface. On the other hand, inspired by recent results on the "universality" of the nodal Fermi velocity 
f i?< (at low energies) in certain cuprates, a prominent role of this quantity is also expected. 

B. The coupling parameters and thermodynamic properties 

Here, we recall that relevant dynamical information is contained in the analytic continuation G{k, cj + zO^) to points 
just above the real frequency axis, known as the "retarded" Green's function [331 . One is therefore led to continue 
the electronic self-energy I](fc,ia;„) analytically by I](fc,a;-|-iO"'") = Ei(fc,a;) -|-iE2(A:, w), where the bare electron band 
energy is determined by the poles of the Green's function G{k,oj + tO+) or the zeros of G~^{k,uj + iO+). Assuming 
that a pole occurs near a; = 0, one gets 



G-\k,uj + iO+) =uj-Sk -Ei(fc,u;) -iE2{k,uj) 

dT.i{k,Lo 



CJ 1 - 



du) 



-[efc + Si(fc,0)]-zI]2(fc,w). 

(4) 



Then, the pole of G occurs at a frequency Wq given by ujq = Ej^ — IJIt}^^ with the quasiparticle scattering time defined 
by T^^ = — 2 (1 — d^Y^x) T,2{k,Ek), and the electron dressed band energy Ek by 

^fc = (i-a^Si)-M£A,- + Si(fc,o)] . (5) 

Now, as a manifestation of the electron-phonon coupling interaction one can introduce the mass renormalisation 
of the electronic dispersion at the energy scale associated with the phonons. This may be technically defined by a 
mass-enhancement parameter A* [401, i.e., 

A^ ^ -d^^iUo ■ (6) 
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FIG. 1: Electron-phonon spectral density a^F{u) for (a) La2-^Sr^Cu04, (LSCO), (b) Bi2Sr2CaCu20s+^ (BSCCO), and (c) 
YBa2Cu30e+x (YBCO), determined by different methods. The sohd hnes correspond to the method in Ref. [42|, dashed lines 
to the method in Ref. [43|, and the dotted line in (b) corresponds to the method of Ref. [43|. In (d) we shown a schematic 
picture of the Migdal-Eliashberg approach from the lowest order Feynman diagram for the electron-phonon interaction. The 
wavy line represents the phonon Green's function D, the double solid line is the renormalised electron Green's function G, and 
the shaded circles represent the small electron-phonon vertex corrections neglected in the ME approach. 



On the other hand, the strength of the electron-phonon interaction is commonly estimated from the so-called boson 
coupling parameter A. This dimcnsionlcss parameter is commonly defined in terms of the electron-phonon spectral 
density as 



A = 2 



^d.^^^m 



(7) 



and it will be eventually related to the superconducting transition temperature. 

We want to emphasise that this quantity is not to be straightforwardly identified with the mass-enhancement 
parameter A* and a necessary distinction between them is essential for the overall description of the available ARPES 
data [34|. In fact, as it was argued in that paper, equality would just be warranted at low temperatures and within 
the Migdal-Eliashberg approximation for the phonon coupling. 

As one can easily appreciate, the spectral density function a^-F(z^) plays a fundamental role in the electron-phonon 
coupling. Since the electron-phonon interaction manifests itself not only in the electronic properties but also in 
phononic properties via the density of phonon states F{i/), the investigation of phonon spectra may provide valuable 
information on the strength and character of electron-phonon coupling. The information about the phonon spectrum 
can be taken from the inversion of the Eliashberg equations in an analysis of tunnelling data |37l . |4l| or from inelastic 
neutron scattering data for superconductors [43, 13 • Along these lines, we have evaluated the electron spectral 
densities for LSCO, BSCCO, and YBCO in the fiat model of Ref. [43 and also by solving the isotropic Eliashberg 
equations on the Matsubara frequencies. Strictly speaking, the structure of these spectral densities is restricted to the 
isotropic nodal direction under the assumption that a^F{i') = G{v) x C, with C an adjustable constant and G{v) the 
generalised phonon density of states extracted from the inelastic neutron scattering experiments [4J| . Our results are 
shown in Fig. 1 (Shiina). In addition, other reproducible methods to calculate the electron-phonon spectral density 
have been taken into consideration, avoiding the intrinsic complexity in the evaluation of the matrix elements that 
determine the spectral density a'^F(v) from ab initio calculations . Sp ecifically, we mean the simple method by Islam 
& Islam 143 for all samples, and the method by Gonnelli et al. [41| in the par ticular case of BSCCO. The method 
of Refs. [42| & [43] is based on the INS experimental data by Renker et al. (44l - l46l |. On the other hand the method 
within Ref. I4l| is based on the tunnelling data reported in the Refs. (47l.l48|. 
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FIG. 2: Plots of the critical temperatures Tc (top), the ratio gap 2Ao/fcsTc (middle) and the gap Ao (bottom) for LSCO, 
BSCCO and YBCO, all of them represented as functions of the electron-phonon coupling parameter A, and correspondin g to 
the spectral densities of Shiina & Nakamura [43[ and Islam & Islam [421 • ^^i addition the spectral density of Gonnelli et al. l4p 
has been considered in the BSGCO case (Fig. [ij. We have used three different approximations, i.e.: McMillan's formula [49| . 
the Allen- Dynes formula [50], and Kresin's formula [5ll |. 



Below, we consider the possibility that the electron-phonon spectral densities presented above (Fig. [T]) allow to 
explain the high- Tc values, and the zero temperature gap observed in experiments. 

In Fig. [2] wc show our results for Tc, the ratio gap 2Ao/fcBTc and the zero temperature gap Ag, all of them based 
on three different approximations: 
(1) The celebrated McMillan's equation ^, 



Tc = Y^ea;p 



-1.04 



1 + A 



A-Ai*(l + 0.62A) 



(8) 



with the phonon characteristic energy defined by uji = 2S/X, S being the area under the spectral density a^F{y), and 

jji* the so-called Coulomb pseudopotential 

(2) The Allen and Dynes formula [SO], which is obtained by replacing loi in Eq. ^ by 



'^log 



A ./n v 







(9) 



In the systems analysed here, we get u>^^f^^^ ~ l6meV, wf"^*-^*^ 



25meV, uj^^^^'^'^ ^ 33meV, wP'^o ~ 39meV, 



(3) Finally, the less conventional Krcsin's formula [51| 



2 



-1/2 



Tc = 0.25 w exp 1 , (10) 

\^eff J 

where tu = [(2/A) ^^ va^F{v)dvYl^ and A^// = (A - m*)[1 + 2^** + (3/2)A/i*ej:p(-0.28A)]. 

As recalled before [SJ], a satisfactory perturbation theory for the Coulomb interaction does not exist. Then, a 
phenomenological point of view is usually adopted, considering that the superconducting components of the Coulomb 
self-energy turn out to have only a small effect over the bare energy of the quasiparticles. In this sense, in a first 
instance, we considered a typical value of the Coulomb's pseudopotential, /i* — 0.13 [34|. In this paper, the Coulomb 
effects and their influence on the electron-phonon interaction picture will be thoroughly discussed. 

As regards the ratio gap, this quantity has been calculated by using the common formula, 2Ao/A:bT'c = 3.53 [1 + 
12.5 (r,/c^iog)2 ln(L^iog/2Tc)]. 

From Fig. [2] several conclusions can be drawn about the relevance of the phonon- mechanism [42|, |43 . Apparently, 
from the theoretical point of view the high values of the critical temperatures strongly depend on the approximation 
used to calculate them, and even in some cases on the inversion method used for obtain the spectral density. In this 
sense, the electron-phonon coupling parameter A has been usually considered a fit parameter used for obtaining the 
correct experimental critical temperatures, and then used to conclude strong or weak infiuence of the electron-phonon 
interaction on the superconductivity properties. Nevertheless, the strong dependence of the derived parameters on 
the particular formulae considered, does not favour this procedure. Then, our position here is determine the electron- 
phonon coupling parameter from experimental available spectroscopies as the energy dispersion curves obtained in 
ARPES experiments and compare it with the Tc value calculated with different approximations an different spectral 
functions. Outstandingly, the clearest point is to check if, regardless the approach used for obtaining the spectral 
function, it leads or not to correct predictions for the observed values of Tc and Aq. 

We should comment that other choices of the spectral density in which the interaction mechanism has a magnetic 
origin are also possible. Nevertheless, in this picture, one must accept that the magnetic mode can be described by the 
same Feynman diagram as the phonon mode (Fig. [TJi). This hypothesis was first proposed by Berk and Schrieffer |52l | 
and recently invoked by Dahm et al. l2p . To support our "phononic choice" we want to recall that the absence of the 
magnetic-resonance mode in LSCO [III, BSCCO over-doped (x=0.23) [53| and its appearance only below Tc in some 
cuprates (e.g., Bi2201 [3y|) are not consistent with the idea that the kink dispersion has a magnetic origin. On the 
other hand, recent studies on electron doped systems [2j, [2a, [2^, [23] have shown that the intensity of the magnetic 
resonance mode is seemingly weak in comparison with the phonon mode to be considered as the cause of the strong 
electron energy dispersion measured by ARPES. However, this issue will be reconsidered in our final discussion. 

III. NODAL ARPES SPECTRA 

It is well known that angle-resolved photoemission spectroscopy (ARPES) has become a powerful technique for 
evaluating the quasiparticle self-energies I](k,aj) in condensed matter, despite the complexities in understanding the 
photoemission processes involved (refer to our discussion in Ref. [34[ and references therein) . On the other hand, within 
the electron-phonon coupling scenario and starting from Eq.(I2]) within the ME approach, the electronic self-energy 
can be obtained from the complex expression 



S(k, iu) = / diya^F{k, v) | -2t 
with the function T{i', uj,T) defined by 



N{v) + \ 



T(i.,w,r) , (11) 



T,,.,r,.*(l..^)-.(1-.^). ,., 

Here ^/'(z) arc the so-called digamma functions. Now, recall that the bare band energy Sk relates to the dressed 
band energy £'k by E\^ = ek -l- Re'E{Ek). Whereas the direct extraction of the self-energy from experiments appears 
to be troublesome because the underlying band structure of the bare electrons is a priori unknown, a theoretical 
determination of the bare band structure and its relation to the full energy renormalisation effects observed in the 
experiments seems much more attractive. In this sense, the nodal ARPES spectra discussed below are of great 
importance to check the validity of quasiparticle concept and understanding the nature of the interactions involved. 
As a description of the ARPES technique and related photoemission spectroscopies is beyond the scope of this 
paper, we consider a pertinent suggestion for the readers the review by A. Damascelli, Z.-X. Shen and Z. Hussain 



[Rcf. [5J]], and the works contained in the most recent books of Stefan Hiifner [Ref. [2|], X. J. Zhou et al. [Ref. y] 
and by J. R. Schrieffer and J. S. Brooks [Ref. ^]. 

The purpose of this section is to go through the prominent resuhs obtained from ARPES measurements on the 
(0,0)-(7r, tt) direction in the Brillouin zone (nodal direction) avoiding the coniphcations of the superconducting gap. We 
want to emphasize that in the nodal direction the rf-wave superconducting gap in cuprates is zero and the interactions 
involving a low-energy excitation appear as a well-defined slope change in the electronic energy-momentum dispersion 
in a similar energy scale (E — Ep ^ 40 — SOmeF or kink \n. la. [lll - [l3l [l8| (see Fig. [3]). On the other hand, the antinodal 
direction denotes the (tt, 0) region where the anisotropic character of the d-wave superconducting gap seems to be 
unavoidable flO. [l3l - [l6{ . In fact, any theoretical assumption to include the d-wave character of the gap function could 
be clouding the real influence of the phonon-coupling modes or even any other mechanism, strong or weak, it involved 
in the superconducting pairing properties. In this sense, is not queer that in the antinodal direction the kink-topology 
is not universal in comparison with the called universal nodal Fermi velocity [8|. Special attention will be payed on 
the electron-phonon interaction in HTSC-cuprate compounds along the lines of the previous discussion. 

Let us to mention that the energy distribution- and momentum distribution-curves (EDC - MDC) are the two most 
popular ways for analysing photoemission data. The dichotomy between the MDC- and EDC-derived bands from the 
same data raises critical questions about its origin and also about which one represents the intrinsic band structure. In 
this sense, we want to recall that at the larger bandwidth along the nodal direction, the MDC method can be reliably 
used to extract high quality data of dispersion in searching for fine structure. It has also been shown theoretically that 
this approach is reasonable in spite of the momentum-dependent coupling if we are only interested in identifying the 
mode energies Q. In a typical Fermi liquid picture, the MDC- and EDC-derived dispersions are identical. Moreover, 
in an electron-boson coupling system the lower and higher energy regions of the MDC- and EDC-derived dispersions 
are still consistent, except right over the kink region [1^. Below, the good agreement between our theoretical data 
and the experimental facts supports this general picture. 

Considering that the kink effect is a common feature of the HTSC cuprates, i.e., it is observable regardless of the 
doping level or the temperature at which the measurement is performed, the only possible scenario seems to be the 
coupling of quasiparticle with phonons. The key issue seems to be the existence of an energy scale (in the range of 
40-80 meV). As in previous work, our position in this paper has been to determine the evolution of the bare electron 
band energy Sk through the systematic evaluation of the dressed band energy Ek as a function of the electron-phonon 
coupling A in several HTSC cuprates, i.e., LSCO, BSCCO and YBCO (Fig.|31). To be specific, we propose to generalise 
use of the equation Ek — Sk + Re{Yi{Ek)} from which one can consider that e^ implicitly depends on Ek through the 
boson coupling parameter A. We recall that £k is not directly available from the experiments. Instead, the electron 
momentum dispersion curve Ek{k — kp) may be measured from the ARPES experiments. This fact, along with an 
ansatz for the ARPES "bare" dispersion will allow to obtain A as a unique constrained parameter that better fits the 
observed kink topology in a particular sample [2l|,|25|. Thus, based on an interpolation scheme between the numerical 
behaviour of the dressed energy band and the experimental data, one can introduce a universal dispersion relation [3J| 
which allows to reproduce in a quite general scheme the nodal dispersions close to the Fermi level in LSCO (Fig. |4]) , 
BSCCO (Fig.©, and YBCO (Fig. [T]) samples. 



As discussed in Ref. [34| . to a first approach, our relation is defined by 



k-kp^-^il-SX) , (13) 

with vp^ the Fermi velocity at low-energies. Recall that, contrary to the behaviour at high energy values, this limit 
of the Fermi velocity is rather independent of the chemical structure and doping levels and can be obtained as the 
slope of the lower part of the momentum dispersion curve Q . Note that then, 5 is the only "/ree" parameter required 
for incorporating the specific renormalisation for a given superconductor. Let us emphasise once again that the 
correct mass-enhancement parameter defined in Eq.|6] cannot be directly identified with the electron-phonon coupling 
parameter A, as well as the complete electron self-energy Si is not to be straightforwardly identified with the electron- 
phonon self energy Sep (though it is an important part of it). As a first approximation to this fact, recall that in 
Eq. (1131) the involved quantities can be assumed to be not far from their values at the Fermi level and consequently, 
one is enabled to replace eu by 

Ek - Ei(Sfc) = {k- kp)vp^ (1 - 5\)-^ . (14) 

As it was shown in Ref. [3J], starting from this relation and Eq. ^ one can show that the parameter 5 is to be 
physically defined by (5 = (A*/A)/(l -|- A*); i.e., (1 — 5\)~^ = 1 + A*. Recall that the proof relies on the fact that A* 
is defined from the electronic group velocity Wk = {l/fijdE/dk, which changes by a factor of 1/(1 -|- A*) close to the 
Fermi level [S^. This allows to write Eq. P^ in the alternative way 
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FIG. 3: Top: The renormalised energy Ek as a function of the bare band energy et assuming that the nodal dispersion is 
a consequence of the electron-phonon couphng. AU curves have been obtained using Eq. (|lip and the EPI spectral density 
a'^F{uj) calculated by the more rigorous method of Ref. ^ (Fig.[T]) in LSCO (left), BSCCO (centre) and YBCO (right). Very 
similar results are found under the use of any of the densities shown in the Fig. [l] The curves are labelled according to the 
e-ph coupling parameter "A" and range between the bare dispersion A = 0.0 (straight line) and about A = 4.0 for LSCO (at 
T=20K), A = 2.82 for BSCCO (at T=70K) and A = 2.74 for YBCO (at T=30K). Bottom: Simulation of the corresponding 
renormalised energies Ek as a function of the momentum k -- kp along the nodal direction from our own model fEQ. (|13[) ) with 
Slsco = 0.185, Sbscco = 0.354, and Sybco = 0.365. 



Comparing both equations, to the lowest order, the dimensionless parameter 6 is basically the ratio between the 
defined mass-enhancement and phonon-coupling parameters S w A*/A. Remarkably, the most general Eq. (jisp allows 
to consider higher-orders in the electron-phonon coupling from the physical meaning of 6, or even consider additional 
mechanisms from the physical meaning of A*. 

In this work, we have extended previous studies to the consideration of different spectral densities, and also analysed 
the implications of different schemes for obtaining Tc- This will allow to check the relevance of the phonon mechanism 
and the relevance of Eq. (jl3p for obtaining the coupling parameter A. The main facts that arise from this analysis are 
developed in following subsections, and summarised in figures S] - [HI 
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FIG. 4: The renormalised electron quasiparticle energy dispersion Ek as a function of the momentum k — kp for several samples 
of La2-i:Sri:Cu04, measured along the (0,0)-(7r,7r) nodal direction at a temperature of 20 K. The doping level x ranges between 
0.03 (top left) up to 0.30 (bottom right). The experimental data (symbols) are taken from Ref. [3| and the theoretical curves 
have been obtained using the spectral densities of Islam & Islam [Ref. [42|] (open squares), and Shiina & Nakamura (solid 
lines) mi.] 



Results in LSCO 



In Fig. |3] we show the singular features involved within the bare electron band energy curves [Ek{sk)] and the 
renormalised dressed electron band energies [Ef^{k—kp)]. One can observe a structure in the renormalised quasiparticle 
energy Ek ranging up to roughly 80 meV. In fact, the bands rapidly approach Ep from high-binding energy and 
suddenly bend at 40-80 meV, showing the kink topology described in the ARPES experiments. 

For quantitative purposes, here, we have considered wf< = 2ey-A as related to the experimental data of Refs. [3- 
0, [ll|. In Fig. IHwe show the results found in La2-xSrxCuOi covering the doping range (0 < a; < 0.3). Remarkably, 
in this material, the hole concentration in the Cu02 plane is well controlled by the Sr content x and a small oxygen 
non-stoichiomctry. Within this range, the physical properties span over the insulating, superconducting, and over- 
doped non-superconducting metal behaviour. Superconducting transition temperatures Tc in the interval of 30-40K 
have been observed by Bednorz and Miiller [l| and others [53, [Sg]. On the other hand, the best fit of the whole set 
of experimental data has been obtained for 5 = 0.185 (Fig. [5]). 

In fig. [5] the evolution of the parameter A as a function of the doping level is shown. Taking advantage of 
the widespread availability of experimental information in LSCO, one can fit the data to the simple expression 
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FIG. 5: Evolution of the e-ph coupling parameter A as a function of the dopant content in La2-xSrxCu4 samples. The A-values 
have been obtained from our best fit with the nodal kink dispersion (black squares) shown in fig. [J] Correspondingly, the 
evolution of the area 5* as a function of the dopant content is also shown (right scale) . 



A = 2ujexp{—jx) + 1, within a precision factor around of 95%. a) is the ratio between the phonon characteris- 
tic energies introduced by McMillan [i^l, wi == {2/X) J^ a'^F(iy)di' = {2/X)S, and Allen and Dynes [501, i^iog = 
exp{{2/X)J^\n{iy)[a^F{iy)/iy]diy}, i.e., w=wi/a;iog. We get Wj^f'^'^ ~ 16.U55meV and wf^'^o ~ 25.2627 meK We 
want to clarify that, although this equation can be just considered as a useful relation between the physical and 
chemical properties of LSCO, the shaping of other HTSC by similar expressions cannot be guaranteed. Regarding the 
critical temperatures, it must be emphasised that for LSCO, the A-values obtained from our model [3J] allow to explain 
the observed superconducting temperatures Tc ^ [53, Ha] . This can be done either by using the McMillan formula 
^^ or in a more general context, by solving the Eliashberg equations [43|. Nevertheless, we want to emphasise on 
how these facts have to be interpreted recalling that several discrepancies could be emerge depending on the formula 
used for obtaining Tc- From the previous section we have concluded that the critical temperature Tc should not be 
considered as a fit parameter for adjusting the theory, i.e.: one should not predict A from the approximate Tc formulas 
and then use it for calculating the electron self-energy. Outstandingly, this attractive idea has led to unfortunate 
underestimates of the phonon contribution to the photoemission kink in HTSC [22, [54I . On the other hand, to our 
knowledge, the most suitable way for determining the infiuence of an interaction mechanism in the pair formation for 
HTSC is (i) evaluate the strength of the boson coupling mode from the electron renormalisation effects and then (ii) 
solve the Eliashberg equations for the superconducting Tc- In addition, a semiempirical approach has been introduced 
in this paper by using the celebrated equations for T'c(A) referred above. From such analysis, we conclude that the 
consideration of the electron-phonon interaction in LSCO strongly suggests that the high Tc values can be caused by 
conventional electron-phonon coupling, in agreement with the conclusion of Weber [531 • I^. that case they obtained 
A = 2 — 2.5 in the range 0.2 > x > 0.1 within the framework of the nonorthogonal-tight-binding theory of lattice 
dynamics, based on the energy band results of Mattheiss [53], and remarkably corresponding to the range of the EPI 
coupling parameter A predicted in our model. Moreover, within our model the quasiparticle parameter "A" calcu- 
lated from the gradient of the DFT self-energy is basically to be identified with the mass-enhancement parameter A*. 
Thus, although the density-functional theory gives a correct ground-state energy, the bands do not necessarily fit the 
quasi-particle band structure used to describe low- lying excitations in agreement with Ref . [2lj , and the experimental 
quasiparticle parameter Xexpt matches the electron phonon coupling A involved in the standard Migdal formalism 
for analysing the experiments [SJ]. On the other hand, it must be emphasised that in the case of Ref. [5^, the Tc 
values were also reproduced. Moderate discrepancies between the predictions of our phenomenological model and the 
analysis of Refs. [53] & [SO] rnay be ascribed to some uncertainty in the experimental spectral densities. Remarkably 
the high Tc values observed in LSCO have been obtained from the McMillan equation, or in a most general way by 
solving the Eliashberg equations according to Ref. (43j . 
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FIG. 6: Same as in Fig. |4] but in samples of Bi2Sr2CaCu208+a;, all measured along the (0,0)-(7r,7r) nodal direction. The 
experimental data for the normal state (full diamond) and superconducting state (full circle) both have been taken from the 
Ref. [131 ■ The theoretical curves (lines -diamond, or -circles) have been obtained from our dispersion relation ( Eq. ll3|l according 
to the best fit values for e-ph coupling parameter A using the spectral densities of Islam & Islam [Ref. [i^l] (top), Shiina & 
Nakamura [Ref. [43|] (middle), and Gonnelli et al. [Ref. [43|] (bottom). The different plots correspond to the doping levels of 
BSCCO: under-doped "x=0.12" (left), optimally doped "x=0.16" (centre), and over-doped "x=0.21" (right) respectively. 



B. Results in BSCCO: analysis of ARPES data 



In Fig. |6]we display the results found for the renormalised electron quasiparticle dispersion Ek as a function of the 
momentum k ^ kp for several samples of the double layer Bi2Sr2CaCu20g^x measured along the (0,0)-(7r,7r) nodal 
direction. BSCCO curves were predicted by our interpolation method and different inversion methods to determine 
the electron-phonon spectral densities (see fig. [J), for the under-doped (UD70) "x ~ 0.12. Tc ~ 70K" , optimally 
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FIG. 7: Same as in Fig.|4]but in samples of YBa2Cu306+a:, all measured along the (0,0)-(7r,7r) nodal direction. The experimental 
data have been taken from Ref. [Ig. The theoretical curves have been obtained from our dispersion relation [Eq. (|13|l ] according 
to the best fit values for e-ph coupling parameter A using the spectral densities of Islam & Islam [Ref. [42|] (open squares), and 
Shiina & Nakamura [Ref. [43]] (solid line). The different plots correspond to the doping levels: "x=0. 4" (under-doped), "x=0.6", 
and "x=0.85" (over-doped) respectively. 



doped (0PD90) "x = 0.16, T^ w 90A'", and over-doped (OVD58) "x = 0.21, T^ « 58K" samples. The experimental 
data were taken from the work by Johnson et al. 17 1 with vp^ — 1.6el/-A as a value consistent with the experimental 
results of Refs. 0, [13, [ij, [l3l- In a first approach, the best fit with experimental data has been found for S = 0.354 
within the same approximations done in the previous section, i.e., for determining the critical temperatures we have 
assumed /i* = 0.13. On the other hand, in this work we show that regardless the method used for obtaining the 
electron-phonon spectral density (fig. [T]) the same conclusions are obtained (fig.©. Recall that, similarly to the case 
of LSCO, our analysis fits well the "A" values predicted by other authors and from very different models |12l. l4ll. |43| 
(see table 2 in Ref. pM| ) . However, as regards to the critical temperatures, only for the BSCCO UD70 sample the 
A values obtained can be consistent with the experimental facts, i. e., from the spectral density of Islam & Islam 
the critical temperatures obtained are: T^^'-<=^™(A = 2.15) = 70.0/-^ and TJ^^MiZian (;^ ^ 2.15) = 64.81A'. From the 
spectral density by Gonnelh et al.: Tf'''"''"{\ = 2.15) = 67MK and T^i<=Miiian^^ ^ 2.15) = 58.19A'. However, we 
must mention that a strong reduction of Tc {'^ 40%) is found from the method by Shiina & Nakamura [43| and the 
Allen & Dynes formula [501 • Perhaps, this fact could be considered as a first signal about the need of considering 
additional perturbation mechanisms previously not involved in the complete solution of the matrix elements for the 
Eliashberg equations, or perhaps and even too, this fact reveals that the flat model used in the Ref. [4j| is not 
consistent with this experimental facts. This will be further discussed in a forthcoming section (Sec lIIID]) . 



C. Results in YBCO: analysis of ARPES data 



Regarding the YBCO family, that is characterised by double Cu02 planes sandwiched between single CuO chains 
a similar analysis has been done, but remarkable difficulties have been found. 

In Fig. [7] we display the results found for renormalised electron quasiparticle dispersion E^ as a function of the 
momentum k — kp for several samples oiYBa2Cu30e+x measured along the (0,0)-(7r,7r) nodal direction. The data were 
collected from our interpolation method and different inversion methods to determine the electron-phonon spectral 
densities (see figH]), in the under-doped "a; = 0.4, T^ « 35A"'(UD35), "x = 0.6, T^ w 61A'"(UD61), and over-doped 
"a; = 0.85, Tc ~ OOAT" (OVD90) samples. The experimental data were taken from the work by Borisenko et al. [I3 with 
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FIG. 8: Plots of the critical temperatures Tc (top), the ratio gap 2Ao/kBTc (middle) and the gap Ao (bottom) for BSCCO 
(top) and YBCO (bottom), all of them represented as functions of the electron- phonon cou plin g parameter A and the Coulomb 
pseudopotential /i*. A has been determined from the spectral densities of Islam & Islam [42| (left), Shiina & Nakamura [43| 
(centre), and Gonnelli et al. [4l[ (right) in the specific case of BSCCO (fig. [T|. We have used three different approximations, 
the Kresin's formula [5l||, the McMillan's formula [40], and the Allen-Dynes formula [SOI ]. 



Vf^ = 1.63ey-A as a value consistent with the experimental results reported by those authors. In the first approach, 
the best fit with experimental data has been found for 6 = 0.365 [3J]. In this work we show that independently of the 
method used for obtaining the electron-phonon spectral density (Fig. [T]) the same conclusions are obtained (fig. [7]). 
Recall that, similarly to the LSCO and BSCCO cases, our analysis fits well the "A" values predicted by other authors 
and from very different models [4j,[53,l60| (see table 3 in [3J]). On the other hand, as regards Tc, we have previously 
concluded that the electron-phonon contribution is very weak considering the low critical temperatures obtained by 
different approaches notwithstanding, the kink effect was entirely reproduced |34{ . A discussion on the mismatch 
between the electron-phonon coupling parameters and the superconducting properties in YBCO is put forward in the 
following section. 
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D. Electron-phonon coupling and superconductivity in BSCCO and YBCO 

The analysis of the results for the BSCCO OPD90 and BSCCO OVD85 cases, already revealed in previous results 
(Ref. (SJI) that the influence of the electron-phonon mechanism is seemingly weak in spite of the fact that the kink 
effect is reproduced within this scenario. In fact, independently of the method used for obtaining the electron- 
phonon spectral density (Fig. [ij the maximal values for the electron-phonon coupling parameter which reproduce 
the kink structure in BSCCO OPD90 (A ~ 1.3) and OVD58 (A ~ 0.85) samples (Fig. i]), both underestimate the 
experimental critical temperatures in around of 56% and 67% respectively. This value has been obtained by assuming 
a conventional value for the Coulomb pseudopotential (/.i* = 0.13) (Fig. [2]). Now, recalling the impossibility of 
describing a satisfactory perturbation theory for the Coulomb effects, here we have analysed the influence of fj,* 
on the Tc approximations to determine possible enhancements of the thermodynamic properties as functions of the 
electron-phonon coupling parameter A and their corresponding spectral densities (Fig. |8]). In particular, taking the 
value fj,* = 0.001 instead of assuming the conventional fj,* = 0.13, the Coulomb effects arc even more negligible 
in the renormalised energy. Thus, in BSCCO OPD90 one obtains an enhancement of around 30% {Tc ~ 68K) as 
compared to the conventional Tc value (Tc ~ 40-fC). This fact is even more notorious in BSCCO OVD58 where Tc 
improves around 37% {Tc '^ AlK). However, even when favoured, the phonon mechanism does not seem to be the only 
mechanism to be considered in the pair formation that leads at the higher critical temperatures observed in BSCCO 
OPD90 and OVD58 samples. Notice that, in order to show this, we have assumed the most favourable scenario 
for the phonon hypothesis, i.e.: 1. we have considered the methods of Islam & islam (Ref. [i^) and Gonnclli et al. 
(Ref. [41|) to determine electron-phonon spectral densities. 2. We have assumed a very weak Coulomb pseudopotential 
(/i* = 0.001), and 3. The empirical formula of Krcsin to determine Tc has been used (see figm-es [T] & [5]) . Within 
this scenario, we get Xisiam{T^'''''''' = 91K) ~ 1.82 and \Gonnem{T^'^'''''' = 91K) ~ 1.93, with the rcnormalisation 
parameters for the ARPES "bare" dispersion S = 0.264 and 6 = 0.248 respectively for BSCCO OPD90. Regarding 
the case of BSCCO OVD58, the electron-phonon coupling parameters obtained are Xisiam{T^^'^^^^ = b^K) ~ 1.12 
and \GonneUi{T^^'^'^™ = 58K) ~ 1.17, with the same S values. Then, we observe a widening of the kink effect 
suggesting the existence of at least one additional perturbation mechanism reducing the momentum of the dispersed 
quasiparticles (see Fig. [51). 

There are several possibilities as candidates of additional perturbation mechanisms, a suitable one being to 
investigate hidden factors in the intensity of the magnetic resonance mode and its evolution with the doping 

level [ni, [n, m, m, [n in . 

Now, we will consider the same hypotheses dealed above for the BSCCO family, i.e.: improving the phonon mech- 
anism by assuming a very weak contribution of the Coulomb pseudopotential {^* = 0.001), in an effort to obtain 
appropriate values of Tc in the case YBCO UD35, with the same value of the electron-phonon coupling A = 0.80. 
Going into detail, we have used the methods by Islam & Islam [i^ and Shiina & Nakamura [4^ to determine the 
electron-phonon spectral densities (Fig. [T]), and then different approaches for Tc have been evaluated (Fig. [2). Thus, 
with the method from Rci.^ we obtain: T/^-^e"™ ~ 37.69, T^i'^^riian ^ 4Q q4^ ^^^ j^ Mien-Dynes ^ 33 35^ and from 
the more rigorous method of Ref.gl] we obtain: T^^'^'''' ~ 36.86, T^^-^Mtiian ^ 35 qj^^ ^^^^ j^AUen-Dynes ^ 28.71. 
Then, under the assumption that the Coulomb effects don't play any important role on the electron self-energy, one 
can conclude that the phonon mechanism could be the most relevant interaction mechanism in the YBCO UD35 
sample. However, regarding the Tc values of the under-doped family, they are still underestimated by a percentage of 
62% in YBCO UD61 and aroimd 80% in YBCO OVD90. 

On the other side, in Fig. [S] we show the maximal influence of the electron-phonon coupling parameters that 
reproduce Tc on the quasiparticle dispersion energies measured by ARPES. Analogously to the BSCCO cases, here 
we have assumed the most favourable scenario for the phonon hypothesis. We obtain \isiam{T^^'^^''^ = &1K) ~ 1.20, 
and XshUna{T^'''"'''' = &IK) ~ 1.23, with the ARPES "bare" dispersion renormalised at (5 = 0.260 and 5 ^ 0.256 
respectively (Fig. ^. For the over-doped YBCO OVD90, the electron-phonon coupling parameters obtained are 
A/s/am(?;^'^''"'" = 90K) ~ 1.83 and XshUna{Tf '''"'''' = 90A') ~ 1.89, with the ARPES "bare" dispersion renormalised 
a,t S ~ 0.116 and S = 0.120 respectively (Fig. [9]). Then, the more remarkable widening of the kink effect is again 
suggesting the existence of at less one additional perturbation mechanism. In fact, in YBCO the importance of an 
additional mechanism seems to be more significant that in BSCCO, and it could be related to the appearance of a 
second kink in the under-doped phase. 

Very recent efforts have been focused on the high-energy part of the ARPES data (0.2-1.5 eV) revealing additional 
changes in the electron energy dispersion ^^waterfalls^^ at '^ 0.4 eV [1^, [20, [2J, [2^, [6l|. The presence of a kink at 
these high energies immediately raises the question: What type of excitation spectrum is required to produce such 
rcnormalisation effect? One possibility is that correlation effects at high orders strongly renormalises the electron- 
phonon coupling. Moreover, with regard at the photoemission process and the inversion of the data measured in 
YBCO samples one should take special care with the residual 3-dimensionality and the possible effect of the CuO 
chains. Nevertheless, these high energy features remain under intense debate [l^, [20, U% \23i, iMi i63J because it is still 
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FIG. 9: The renomialised energy Ek as a function of the momentum k — kp under the assumption that the nodal dispersion is 
a consequence of the electron-phonon couphng and the Coulomb effects are completely contained in the bare energies. Then, 
we assume the very weak Coulomb pseudopotential /i* = 0.001 within the more favourable approximation for the calculus of Tc 
within t he p honon hypothesis, i.e. the Kresin Formula (Fig. El) and the EPI spectral density oi^ F(uj) obtained by the methods 
of Refs. Ha Ell in BSCCO cases (top), and Refs. ^^ in YBCO cases (bottom). 



unclear whether they represent intrinsic band structure or not [if 



IV. CONCLUDING REMARKS 



In summary, it has been shown that, regarding the ARPES measurements in superconducting cupratcs the strong 
renormalisation of the band structure, customarily related to the dressing of the electron with excitations can be 
explained in terms of the conventional electron-phonon interaction. In fact, our results suggest that the electron- 
phonon interaction strongly influences the electron dynamics of the high- Tc superconductors, and it is an important 
mechanism linked with the Fermi surface topology. Thus, the electron-phonon interaction (strong or weak) must be 
included in any realistic microscopic theory of superconductivity. 

Remarkably, the proposal of a simple linear dispersion relation for the electron bare band energy, i.e., Sk ~ {k ~ 
kF)vF<{^ + A*), with the electronic mass-enhancement parameter A* defined self-consistently by the relation 1 -|- A* = 
(1 — dX)~^ allows to reproduce the appearance of the ubiquitous nodal kink for a wide set of ARPES experiments 
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in the cupratcs. Here, A has been interpreted as the strength of the electron-phonon interaction or electron-phonon 
couphng parameter, and the dimensionless parameter d recognised as a universal property for each family of cuprates. 
In addition, at the lowest order of the Feynman graphs for the electron-phonon interaction this method reassembles 
the "A" values obtained from several discrepant models. 

In this work, we have investigated the influence of the electron-phonon coupling mechanism through different doping 
levels in several families of HTSC-cuprates. On the one hand, we have evaluated different methods for obtaining the 
electron-phonon spectral densities and their influence on the electron bare band energy, and on the other hand, several 
approaches have been recalled to obtain the critical temperatures. Then, our results suggest that at least in the LSCO 
family, and in the so-called BSCCO UD70 and YBCO UD35 superconductors, the electron-phonon interaction could 
be the most relevant mechanism involved in the pair formation that leads to the superconductivity. Our conclusion is 
supported by the experimental evidence of a mass renormalisation of the electronic dispersion curves measured along 
the nodal direction in ARPES and the reported Tc values in good agreement with our theoretical predictions. An 
excellent comparison between the theory and the available collection of experiments is achieved. Furthermore, we have 
evaluated the consequences of assuming an enhanced phonon mechanism, through the reduction of the Coulomb's 
pseudopotential weight. When appropriate Tc values are obtained by this method, a remarkable widening of the 
predicted kink effect arises. This fact, suggest that independently of the approximations invoked and even avoiding 
the influence of the d-wave superconducting gap through the nodal ARPES measurements, it doesn't seem possible to 
elude the existence of additional mechanisms that reduce the momentum of the dispersed quasiparticles in comparison 
with the phonon mechanism. In this sense, despite the fact that in LSCO the influence of the magnetic mode seems 
not relevant, it is not possible to ignore its importance over the electron properties of other HTSC families having in 
mind the recent observation by G. Yu et al. |64l | on a possible connection between the magnetic excitations and the 
superconductor gap for a wide range of materials. 
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